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C3 Abstract 



Preliminary results^] on a measurement of the proton structure function F 2 (x, Q 2 ) 
are reported for momentum transfers squared Q 2 between 1.5 GeV 2 and 5000 GeV 2 
and for Bjorken x between 5 • 10~ 5 and 0.32, using data collected by the HERA 
experiments HI and ZEUS in 1994. F 2 increases significantly with decreasing x, 
qq \ even in the lowest reachable Q 2 region. The data are well described by a Next 

to Leading Order (NLO) QCD fit, and support within the present precision that 
. the rise at low x within this Q 2 range is generated " radiatively" via the DGLAP 

\ evolution equations. Prospects for future structure function measurements at 

^sO ' HERA are briefly mentioned. 
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1 Introduction 



The measurement of the inclusive deep inelastic lepton-proton scattering cross section 
has been of great importance for the understanding of the substructure of the pro- 
ton M. Experiments at HERA extend the previously accessible kinematic range up to 
very large squared momentum transfers, Q 2 > 5 ■ 10 4 GeV 2 , and to very small values of 
Bjorken x < 10~ 4 . The first observations showed a rise of the proton structure function 
F 2 (x,Q 2 ) at low x < 10~ 2 with decreasing x [TJ, [J, based on data collected in 1992. 
This rise was confirmed with the more precise data of 1993 [[§, ||. Such a behaviour 
is qualitatively expected in the double leading log limit of Quantum Chromodynam- 
ics ||. It is, however, not clarified whether the linear QCD evolution equations, as the 
conventional DGLAP evolution in InQ 2 and/or the BFKL evolution || in ln(l/x), 
describe the rise of F 2 or whether there is a significant effect due to nonlinear parton 
recombination 0. Furthermore, it is unclear whether this rise will persist at low Q 2 , 
say of the order of a few GeV 2 . For example Regge inspired models expect F 2 to be 



1 Invited talk given at the XXVI" 1 Workshop on High Energy Physics, Gravitation and Field 
Theory, held in Protvino, Russia, June 1995. The results presented here are reflecting the status of 
the measurements reported by the HI and ZEUS collaboration at the 1995 summer conferences. 



flatter for small Q 2 . The quantitative investigation of the quark-gluon interaction 
dynamics at low x is one of the major goals of HERA. It requires high precision for 



2xF l (x, Q 2 ) + 2(1 - y)F 2 (x, Q 2 ) - ( 2y - *- 1 xF 3 (x, Q 2 ) 



1) 



the F 2 measurement and a detailed study of the hadronic final state behaviour JT(J . 

The structure functions F 1 (x, Q 2 ), F 2 (x, Q 2 ) and F 3 (x, Q 2 ) are related to the inclu 
sive lepton-photon cross-section 

d 2 a e+ P _ 2na 2 
dxdQ 2 ~ xQ 4 _ 

and depend on the squared four momentum transfer Q 2 and the scaling variable x. 
These variables are related to the inelasticity parameter y and to the total squared 
centre of mass energy of the collision s as Q 2 = xys with s = 4E e E p . However, at low 
Q 2 xF 3 can be neglected and the previous expression can be rewritten as a function of 
F 2 and R 

R - F *- 2xF i = Fl ( 2 ) 
2xF x 2xFi y ' 

R(x, Q 2 ) could not be measured yet at HERA, but can be computed, supposing that 
perturbative QCD hold, for a given set of parton densities. Thus F 2 (x, Q 2 ) can be de- 
rived from the double differential cross-section d 2 a / 'dxdQ 2 after experimental and QED 
radiative corrections. The structure function xF 3 (x, Q 2 ) has not been measured yet, 
due to lack of statistics at high Q 2 . However, simple differential cross-section da/dQ 2 
both on neutral currents (exchange of a 7 or Z°) or in charged current (exchange of 
H /± ) have already been published WU ITJJ. In the rest of this paper we will consider 
only the case of the 7 exchange. 

In 1994 the incident electron energy was E e = 27.5 GeV and the proton energy 
was E p = 820 GeV. The data were recorded with the HI Jl3| and ZEUS [14| detectors. 



A salient feature of the HERA collider experiments is the possibility of measuring 
not only the scattered electron but also the complete hadronic final state, apart from 
losses near the beam pipe. This means that the kinematic variables x, y and Q 2 
can be determined with complementary methods which experimentally are sensitive 
to different systematic effects. The comparison of the results obtained with different 
methods improves the accuracy of the F 2 measurement. A convenient combination of 
the results ensures maximum coverage of the available kinematic range. 

In this paper after a description of the data samples (section 2) and of the kine- 
matic reconstruction/event selection used (section 3) we provide the F 2 measurement 
in section 4. and its interpretation at low Q 2 and in terms of perturbative QCD in 
section 5, before giving some prospects in conclusion. 



2 Data Samples 

In 1994 both experiments have reduced the minimum Q 2 at which they could measure 
F 2 using several techniques. For DIS events at low Q 2 the electron is scattered under a 



large angle 9 e ( the polar angles 9 are denned w.r.t the proton beam direction, termed 
"forward" region). Therefore the acceptance of electrons in the backward region has 
to be increased or the incident electron energy to be reduced to go down in Q 2 . This 
was realized as follows. 

i) both experiments were able to diminish the region around the backward beam 
pipe in which the electron could not be measured reliably in 93, thus increasing the 
maximum polar angle of the scattered electron (cf jl^, [L6| for details). This large 
statistic sample, taken with the nominal HERA conditions is called the "nominal ver- 
tex" sample. Its integrated luminosity is between 2 and 3 pb -1 , depending on the 
analysis / experiment . 

ii) Following a pilot exercise performed last year, 58 nb _1 of data was collected for 
which the interaction point was shifted by +62 cm, in the forward direction, resulting 
in an increase of the electron acceptance. This sample is refered to as the "shifted 
vertex" data sample. In HI the low Q 2 region was also covered by analyzing events 
which originated from the "early" proton satellite bunch, present during all periods of 
the HERA operation, which collided with an electron bunch at a position shifted by 
+63 cm. These data, refered to as the "satellite" data sample amount to ~ 3% of the 
total data and correspond to a total "luminosity" of 68 nb _1 . 

iii) Both experiments used DIS events which underwent initial state photon radi- 
ation detected in an appropriate photon tagger to measure F 2 at lower Q 2 (so called 
"radiative" sample |T7|]). The incident electron energy which participate in the hard 
scattering is thus reduced, and so is the Q 2 . 

The luminosity was determined from the measured cross section of the Bethe Heitler 
reaction e~p — > e~pj, measuring the hard photon bremsstrahlung data only. The 
precision of the luminosity for the nominal vertex position data amounts to 1.5%, i.e. 
an improvement of a factor 3 w.r.t the analysis of the 1993 data. For the shifted vertex 
data the luminosity uncertainty is higher (4.5% for HI). The luminosity of the satellite 
data sample was obtained from the measured luminosity for the shifted vertex data 
multiplied by the efficiency corrected event ratio in a kinematic region common to both 
data sets. The uncertainty of that luminosity determination was estimated to be 5%. 



3 Kinematics and Event Selection 

The kinematic variables of the inclusive scattering process ep — > eX can be recon- 
structed in different ways using measured quantities from the hadronic final state and 
from the scattered electron. The choice of the reconstruction method for Q 2 and y 
determines the size of systematic errors, acceptance and radiative corrections. The 



basic formulae for Q 2 and y used in the different methods are summarized below, x 
being obtained from Q 2 = xys. For the electron method 



hj e L 1 — y e 

The resolution in Q 2 is 4% while the y e measurement degrades as l/y e and cannot be 
used for y e < 0.05. In the low y region it is, however, possible to use the hadronic 
methods for which it is convenient to define the following variables 

s = y Z E h-Pz, h . P h T 2 = (E^) 2 + (EivO 2 (4) 

h h h 

Here E,p x ,p y ,p z are the four-momentum vector components of each particle and the 
summation is over all hadronic final state particles. The standard definitions for y h 
and 9h are 

Vh= 2E e tan Y = pJ (5) 

The combination of yh and Q 2 defines the mixed method which is well suited for medium 
and low y measurements. The same is true for the double-angle method which makes 
use only of 9 e and 9 h and is thus insensitive to the absolute energy calibration: 

tan% „ „ cot% 

Vda = a — <r Qda = 4 ^ e 3 s- 6 

tan | + tan^ DA e tanf + tan ^ 

The formulae for the E method are constructed requiring Q 2 and y to be indepen- 
dent of the incident electron energy. Using the conservation of the total E — P z = 
J2i Ei — p Z;i , the sum extending over all particles of the event, 2E e is replaced by 
S + E' e (l — cos^ e ) which gives 

^ ~ £ + ^(l-cos£ e ) Qs " (7) 

By construction y^, and Q| are independent of initial state photon radiation. With 
respect to yu the modified quantity ys is less sensitive to the hadronic measurement at 
high y, since the E term dominates the total E — P z of the event. At low y, yh and y^, 
are equivalent. 

HI measures F2 with the electron and the E method and after a complete consis- 
tency check uses the electron method for y > 0.15 and the E method for y < 0.15. 
ZEUS measures F 2 with the electron and the double angle method and for the final 
results uses the electron method at low Q 2 and the double angle method for high Q 2 
(Q 2 > 15 GeV 2 ). The binning in x and Q 2 was chosen to match the resolution in 
these variables. It was set at 5(8) bins per order in magnitude in x(Q 2 ) for the HI 



experiment and about twice as coarse for the ZEUS experiment in x (similar to HI in 

Q 2 ). 



The event selection is similar in the two experiments. Events are filtered on-line us- 
ing calorimetric triggers which requests an electromagnetic cluster of at least 5 GeV not 
vetoed by a trigger element signing a beam background event. Offline, further electron 
identification criteria are applied (track-cluster link, shower shape and radius) and a 
minimum energy of 8(11) GeV is requested in ZEUS(Hl). HI requests a reconstructed 
vertex within 3a of the expected interaction position, while ZEUS requires that the 
quantity 5 = £ + E' e (l — cos 9) satisfies 35 GeV < 5 < 65 GeV. If no particle escapes 
detection, 5 = 2E = 55 GeV, so the 5 cut reduces the photoproduction background 
and the size of the radiative corrections. The only significant background left after the 
selection comes from photoproduction in which a hadronic shower fakes an electron. 
In HI for instance, It has been estimated consistently both from the data and from 
Monte Carlo simulation and amounts to less than 3% except in a few bins where it can 
reach values up to 15%. It is subtracted statistically bin by bin and an error of 30% is 
assigned to it. 

The acceptance and the response of the detector has been studied and understood 
in great detail by the two experiments: More than two millions Monte Carlo DIS events 



were generated using DJANGO [18] and different quark distribution parametrizations, 



corresponding to an integrated luminosity of approximately 20 pb 1 . The program 



is based on HERACLES [jnj for the electroweak interaction and on LEPTO [20] to 
simulate the hadronic final state. HERACLES includes first order radiative correc- 
tions, the simulation of real bremsstrahlung photons and the longitudinal structure 



function. The acceptance corrections were performed using the GRV [34| or the MRS 



parametrization [[Hp, which both describe rather well the HERA F 2 results of 1993 
for Q 2 > 10 GeV 2 . LEPTO uses the colour dipole model (CDM) as implemented in 



ARIADNE [22 1 which is in good agreement with data on the energy flow and other 



characteristics of the final state as measured by HI p3| and ZEUS pij. For the es- 



timation of systematic errors connected with the topology of the hadronic final state, 



the HERWIG model |25| was used in a dedicated analysis. Based on the GEANT 
program [26] the detector response was simulated in detail. After this step the Monte 
Carlo events were subject to the same reconstruction and analysis chain as the real 
data. 



4 Structure Function Measurement 

The structure function F 2 (x,Q 2 ) was derived after radiative corrections from the one- 
photon exchange cross section 



d 2 a 2na 2 y 



2 



(2-2y + -^—)F 2 (x,Q 2 ) (8) 



dxdQ 2 Q 4 x ° 1 + R 



Effects due to Z boson exchange are smaller than 2%. The structure function ratio 
R = F 2 j1xF\ — 1 has not been measured yet at HERA and was calculated using 
the QCD relation JZ7| and the GRV structure function parametrization. Compared 
to the 1993 data analyses 0, |5| the F 2 measurement has been extended to lower Q 2 
(from 4.5 GeV 2 to 1.5 GeV 2 ), and to lower and higher x (from 1.8 • 10~ 4 — 0.13 to 
5 • 10~ 5 — 0.32). The determination of the structure function requires the measured 
event numbers to be converted to the bin averaged cross section based on the Monte 
Carlo acceptance calculation. All detector efficiencies were determined from the data 
using the redundancy of the apparatus. Apart from very small extra corrections all 
efficiencies are correctly reproduced by the Monte Carlo simulation. The bin averaged 
cross section was corrected for higher order QED radiative contributions and a bin 
size correction was performed. This determined the one-photon exchange cross section 
which according to eq.|| led to the values for F 2 (x, Q 2 ). 

Due to the different data sets available: "nominal vertex" data, "radiative events", 
"shifted vertex" data and "satellite" data, which have different acceptances and use 
for a given Q 2 ,x point different parts of the detectors, cross checks could be made in 
kinematic regions of overlap. The results were found to be in very good agreement 
with each other for all kinematic reconstruction methods used. 

The large available statistics allows to make very detailed studies on the detector 
response: efficiencies and calibration. As a result the systematic errors on many ef- 
fects are reduced, compared to the 1993 data analysis. Here only a brief summary of 
these preliminary errors is given, refering the reader to the original and forthcoming 
F 2 publications: For the electron method the main source of error are the energy cal- 
ibration (known at 1.5% level in 1994), the knowledge of the electron identification 
efficiency and to a lesser extent the error on the polar angle of the scattered electron, 
(59 =lmrad) in particular at the lowest Q 2 and the radiative corrections at low x. For 
the E method, the knowledge of the absolute energy scale for the hadrons, the fraction 
of hadrons which stay undetected in particular at low x, due to calorimetric thresholds 
and to a lesser extent the electron energy calibration are the dominating factors. For 
the double angle method, the major problem comes from the precision in the resolu- 
tion of the hadronic angle at low x and low Q 2 . Further uncertainties common to all 
methods (selection, structure function dependance etc.) were also studied. The pre- 
liminary error on the 1994 data ranges between 10 and 20% with expected final values 
for publication below 10%. 

The preliminary results of HI and ZEUS are shown in fig. [l|, in bins of fixed Q 2 . The 
rise of F 2 at low x is confirmed with the higher precision, and is now observed down 
to the lowest Q 2 measured (1.5 GeV 2 ). The observed good agreement between HI and 
ZEUS and the smooth transition between the HERA and the fixed target (E665, NMC) 
data consolidates this result which can thus be confronted to theoretical expectations. 
The steepness of the low x rise increases visibly with the Q 2 , a characteristic expected 
from perturbative QCD. This rise cannot be attributed to the presence of " diffract ive" 
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Figure 1: Preliminary measurement of the proton structure function F 2 (x,Q 2 ) as 
function of x in different bins of Q 2 . The inner error bar is the statistical error. The 
full error represents the statistical and systematic errors added in quadrature. 



events in the DIS sample since their proportion has been shown to stay essentially 
constant (10%) independently of x and Q 2 |28], |29 . 



5 Low Q 2 and Perturbative QCD 



A test of perturbative QCD is displayed in fig. |2| which represents the results of Next to 
Leading Order (NLO) QCD fit performed by HI as explained in []30| on the data with 
Q 2 < 5 GeV 2 . In order to constrain the structure function F 2 at high x, data from 
the fixed target scattering experiment NMC [EHI] and BCDMS |3z 



are used, avoiding 

regions where higher twist and target mass effects could become important. 
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Figure 2: Preliminary measurement of F 2 (x, Q 2 ). The HI data is consistent with the 
fixed target experiments BCDMS and NMC. The curve represents the NLO QCD-fit 
decribed in the text. The gap visible around 100 GeV 2 corresponds to a boundary 
region between two calorimeters in the HI detector, which is not completely analyzed 
yet. 

The F 2 behaviour can be well described by the DGLAP evolution equations within 
the present preliminary errors. The data for Q 2 values below 5 GeV 2 are also compatible 
with the extrapolation of the fit in this region, as can be seen in the figure. 



This preliminary result is consistent with the published QCD analysis of HI and 
ZEUS 1993 data, which allowed to determine the gluon density in the proton, and 
observe its steep rise at low x as displayed in fig. |3| for Q 2 = 20 GeV 2 . 
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Figure 3: Preliminary measurement of the gluon density from NLO QCD fit. The HI 
result is shown with the complet systematic error band, for ZEUS only the value of the 
gluon density is shown, the error band is similar to the HI one. 

Another test of perturbative QCD lies in observing the asymptotic behaviour as 
suggested by early studies |J . Ball and Forte |33j have recently shown that evolving a 
flat input distribution at some Q , of the order of 1 GeV 2 , with the DGLAP equations 
leads to a strong rise of F 2 at low x in the region measured by HERA. An interesting 
feature is that if QCD evolution is the underlying dynamics of the rise, perturbative 
QCD predicts that at large Q 2 and small x the structure function exhibits double 
scaling in the two variables a and p defined as: 



a 



y/\0g(x Q /x) ■ log(t/t ), P 



logpo/x) 
\J log(iAo) 



with t = log(<2 2 /A 2 



(9) 



In figure § the HI data are presented in the variables a and p, taking the boundary 

- ' - ■ ~" - n - • (4) - 



conditions to be xq = 0.1 and Ql = 0.5 GeV 2 , and A^q = 185 MeV. In a previous 



analysis |Kj the value Ql = 1 GeV 2 was chosen, but the new low Q 2 data seems to 



indicate that F 2 is not yet flat for this Q 2 value. The measured values of F 2 are rescaled 
by 

R' F (a,p) = 8.1 exp [5- + \ log(a) + log^)) , (10) 
\ P 2 1 ) 

to remove the part of the leading subasymptotic behaviour which can be calculated in 
a model independent way; log(i?^F 2 ) is then predicted to rise linearly with a. Scaling 



in p can be shown by multiplying F 2 by the factor R F = R' F e~ 2 " fa '. Here 7 = 2y3/6o 
with b Q being the leading order coefficient of the (3 function of the QCD renormalization 
group equation for four flavours, 5 = 1.36 for four flavours and three colours. 
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Figure 4: Tie rescaled structure functions RfF 2 and R' F F 2 (preliminary) plotted 
versus the variables p and a defined in the text. Only data with p > 1.5 are shown in 
b. 

Fig. |a shows RfF 2 versus p. Scaling roughly sets in for p > 1.5. Fig. [|b, for 
p > 1.5, shows scaling behaviour, namely a linear of log(R' F F 2 ) with o. 

These observations suggest that perturbative QCD could be already valid at Q 2 = 
1 or 2 GeV 2 . Indeed within the present precision we can observe in fig. [| the valid- 



ity of the Gliick, Reya, Vogt (GRV) model |34| which assumes that all low x partons 
are generated "radiatively" starting from a very low initial Q 2 = 0.34 GeV 2 scale, in 
which both gluon and quark densities are "valence" like. This result appears surprising 
since perturbative QCD does not apply at such a low scale, but the HERA results and 
the E665 |35| preliminary results follow the GRV expectations as early as 0.8 GeV 2 . 
More precise data are needed to further constrain the model and draw definite con- 
clusions on the dynamics underlying the low x rise. Nevertheless these results appear 




Figure 5: Preliminary measurement of the proton structure function ^(x, Q 2 ) by HI 
and ZEUS, compared to the results of the E665 and NMC experiments and to the 
"prediction" of the GRV model J3^/ over the full Q 2 range. 



already very promising for the DGLAP evolution equations which might not need to 
be supplemented by the BFKL evolution at low x, in the HERA kinematic domain. 

Focusing now on low Q 2 , the persistent rise of F 2 at low x, when going down in Q 2 
indicates that the photoproduction regime has not been reached yet. This can be seen 
in fig. H which display the behaviour of the total cross-section of the proton-virtual 
photon system as a function of W, the invariant mass of the system (at low x, 
W ~ \JQ 2 /x). F2 is related to the total cross-section of the proton- virtual photon 
interaction ((Jtot(l*p)) y i a 

A 2 

a tot ( 1 * P )^^F 2 (W,Q 2 ). (11) 
The Otot growth can be contrasted with the weak rise with W of the total real photo- 



production cross-section in the same range of W : 20-250 GeV [36, 37 
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Figure 6: Preliminary measurement of the cross-section (J to t{l*P) as a function ofW 
and Q 2 . Results from DIS are compared with the measurements in photoproduction. 
For the readability of the plot, not all Q 2 bins are shown. 



The Regge inspired models DOLA [38 and CKMT|39| which can describe the 
behaviour of (y to t{lP) predicts a rather flat behaviour of F 2 at a few GeV 2 . As shown 
in fig. [7|, the DOLA model clearly fails before 1.5 GeV 2 , while the CKMT model 
which assumes that the "bare" pomeron visible at high Q 2 has a higher trajectory 
intercept (~ 0.24) than the "effective" pomeron involved in "soft" interactions (~ 0.08) 
undershoot the data in a less critical manner. In this same plot we can also notice the 
similarity above 5 GeV 2 between the different parametrizations (GRV, MRS, CTEQ 
f34| , [|0| ) which use essentially the same data to determine their parton distributions 
at the reference scale. 




Figure 7: Preliminary measurement of the proton structure function F 2 (x, Q 2 ) in 
the low Q 2 region by HI and ZEUS, together with results from the E665 experiment. 
Different predictions for F 2 are compared to the data. The DOLA and CKMT curves 
are only shown for the upper row of Q 2 bins; CTEQ3M, MRSG and MRSA ' are shown 
for the lower row; GRV is shown for the full range. 



6 Prospects for Structure function measurements 
at HERA 



The HERA structure function program is still in its infancy, but has already provided 
exciting results at low x. The dynamics underlying the behaviour of the structure 
function can be studied in an exclusive way using jets or particle spectra, since HERA 
is a collider equipped with (two) Att detectors. In the next 2 or 3 years Fl will be 



measured, by taking data at different beam energies in order to keep x and Q 2 constant 
while varying y, thus improving the knowledge on F 2 and on QCD. The statistics will 
increase in such a way that a first measurement of xF 3 will be made, and a precise 
determination of a s should be possible. In 95, both experiments have upgraded their 
detector in the backward area in order to reach lower Q 2 (~ 0.1 GeV 2 ) with good 
precision. This year will thus be devoted to understand the questions raised in this 
paper concerning the low x and low Q 2 dynamics and to open up further stringent test 
of QCD, in particular about the behaviour of the high parton density. 
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